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Abstract 

Josephson junctions with ferromagnetic barrier can have positive or negative critical current 
depending on the thickness dp of the ferromagnetic layer. Accordingly, the Josephson phase in 
the ground state is equal to (a conventional or junction) or to tt (tt junction). When and tt 
segments are joined to form a "0-7T junction" , spontaneous supercurrents around the 0-tt boundary 
can appear. Here we report on the visualization of supercurrents in superconductor-insulator- 
ferromagnet-superconductor (SIFS) junctions by low-temperature scanning electron microscopy 
(LTSEM). We discuss data for rectangular 0, tt, 0-tt, 0-tt-0 and 20 x (0-7T-) junctions, disk-shaped 
junctions where the 0-tt boundary forms a ring, and an annular junction with two 0-tt boundaries. 
Within each or tt segment the critical current density is fairly homogeneous, as indicated both 
by measurements of the magnetic field dependence of the critical current and by LTSEM. The tt 
parts have critical current densities j'J up to 35 A/cm 2 at T = 4.2 K, which is a record value for 
SIFS junctions with a NiCu F-layer so far. We also demonstrate that SIFS technology is capable 
to produce Josephson devices with a unique topology of the 0-tt boundary. 

PACS numbers: 74.50.+r, 85.25.Cp 74.78.Fk 68.37.Hk 
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I. INTRODUCTION 



As predicted more than 30 years ago[T], Josephson junctions can have a phase drop of 7r 
in the ground state. Such it junctions are now intensively investigated, as they have a great 
potential for applications in a broad range of devices ranging from classical digital circuits [21 
El SI E] to quantum bits[6j EJ [9]. Nowadays, 7r Josephson junctions can be fabricated by 
various technologies, including junctions with a ferromagnetic barrier [10] El H21 EH HU [15], 
[16j [T71 [18] , quantum dot junctions[19] [20] [21] and nonequilibrium superconductor-normal 
metal-superconductor Josephson junctions [221 [231 [21] 

In the simplest case the supercurrent density j s across the junctions is given by the first 
Josephson relation 



with the critical current density j c > for a junction and j c < for a it junction. Here, 
is the gauge invariant phase difference of the superconducting wave function across the 
junction (Josephson phase). 

Particularly superconductor-insulator-ferromagnet-superconductor (SIFS) junctions [HI 
Uni HE] are promising since, in contrast to other types of it junctions, they exhibit only small 
damping at low temperatures, which is necessary to study Josephson vortex dynamics as 
well as to use them as active elements in macroscopic quantum circuits. 

Now consider a junction in the x-y plane, which has a region with critical current density 
j® > (0 region) and another region having j£ < (n region). For the sake of simplicity 
let us assume that the boundary between and 7r regions runs along the y direction. When 
is different from or tt the supercurrents flow in opposite directions on the two sides of 
the 0-7T boundary, forming a vortex, with its axis coinciding with the 0-tt boundary (along 
the y direction), that carries a magnetic flux $ = ±$o/2 ($o ~ 2.07 x 10~ 15 Wb is the flux 
quantum) [221 EH1 EZ]- This is true if the junction length L in x direction is much larger 
than the Josephson penetration depth 



Here fiod' is the inductance per square (with respect to in-plane currents) of the supercon- 
ducting electrodes forming the junction. For junctions having electrode thicknesses larger 
than the London penetration depth Al, d! ~ 2A^. Experimentally, such semifluxons have 



j s = j c sin 



(1) 




(2) 



3 



first been studied in the context of cuprate grain boundary junctions|28, 29J or zigzag ramp 
junctions between Nb and YE^CuaOy [SO]- Here, the sign change of the ci-wave order 
parameter of the cuprates leads to the formation of O-ir facets. In junctions with a ferro- 
magnetic barrier the value (and the sign) of the critical current density crucially depends on 
the thickness dp of the F-layer [TU [16]. A junction consisting of various and 7r segments 
can, thus, be formed by selectively etching the F-layer to produce two thicknesses d° F and d F 
of the F-layer such that they correspond to critical current densities j® and with opposite 
signs and j° ~ |jj|[31j. 

In the cuprate/Nb zigzag junctions [321 E21 E31 EI] the facets should be oriented along 
the crystallographic a and b axes of the cuprate electrode, imposing certain topological 
limitations to the 0-tt boundary. In contrast, the SIFS technology allows almost any 2D 
shape of the 0-7r boundary and therefore offers a higher degree of design flexibility. Below, we 
show an example where this boundary forms a loop. Even intersecting 0-7T boundaries should 
be feasible, e.g., by arranging and tt regions in a checkerboard pattern. Unfortunately, 
the present SIFS technology based on a NiCu ferromagnetic layer produces a maximum 
which is much lower than j c ~ lkA/cm 2 of standard Josephson tunnel junctions. Although 

at T = 4.2 K has been increased from some mA/cm 2 for the first junctions [TT], to a few 
A/cm 2 in Ref. [161 and to about 35 A/cm 2 in the present paper, the value of Aj oc 1/y 
is still above 50 /im. Thus, the study of a multi semifluxon system would thus require 
unreasonably large (mm sized) junctions. 

Nonetheless, also (multifacet) junctions with length L < Aj are interesting. For example, 
one can consider an array of many alternating and 7r segments along x, where the lengths of 
individual segments are much smaller than Aj. Such a structure is similar to short multifacet 
cuprate/Nb zigzag junctions[30l [321 [33] or high angle grain boundaries in high T c cuprates 
[35] and can e.g. be used to realize a ip junction — a junction having a phase ±ip in the 
ground state and many other interesting properties [Ml E3 EB1 ESI I4D] . 

The goal of this work is to realize Josephson junctions with various arrangements of 
and 7r segments in order to demonstrate that also complex structures are feasible. We 
characterize these junctions by measurements of current voltage (I—V) characteristics, by 
I C {B) and by low-temperature scanning electron microscopy (LTSEM) |41j. By analyzing 
I C (B), in principle one obtains information on the suercurrent flow and (in)homogeneity 
of the critical current; however, the analysis at least of the more complex SIFS structures 
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may require to consider many unknown parameters (gradients in critical current density, 
local inhomogeneities etc.), making conclusions ambiguous. We thus put a strong focus on 
LTSEM which allows direct imaging of the supercurrent density distribution in the junctions 
(including counterflow areas induced by the 0-7T- segments), close to I c [34J. 

The paper is organized as follows. In Sec. |TT] we discuss the sample fabrication and 
measurement techniques. The experimental results are presented and compared with the 
numerical simulations in Sec. |III[ Different subsections are devoted to various geometries, 
(0 junction for reference, 0-7r and 0-7T-0 junctions, a junction consisting of 0-7r regions pe- 
riodically repeated 20 times, a disk shaped structure where the 0-7T boundary forms a ring 
and an annular junction containing two 0-tc boundaries). All investigated samples are in the 



short limit (L < 4Aj). Finally, Sec. IV concludes this work. 



II. SAMPLES AND MEASUREMENT TECHNIQUES 



A. Sample fabrication 

The Nb|Al 2 03|Ni 06 Cuo.4|Nb heterostructures used for our studies were fabricated, as 
described in Refs. ED H21 In brief, one starts with a Nb|Al 2 03 bilayer (Nb thickness is 
120 nm) as for usual Nb based Josephson tunnel junctions. The thicknesses of the following 
F-layer must be chosen very accurately to realize and ir regions with approximately the 
same critical current density. To achieve that, first the Nio.6Cuo.4 F-layer is sputtered onto 
the wafer with a thickness gradient along the ^-direction to achieve a wedge-like NiCu layer. 
Later on, a set of structures extending along x and consisting of the 0-tt devices to be 
measured, plus purely and n reference junctions, is repeated several times along the y- 
direction. One of the sets will have the most suitable F-layer thickness to yield 7r coupling 
with roughly optimal critical current density. In this way the number of wafer runs which 
are required to get appropriate 0-ir junctions is minimized. After the deposition of a 40 nm 
Nb cap-layer and lift-off one obtains a complete SIFS stack, however without steps in the 
thickness of the F-layer yet. To produce such steps, the parts of the structures that shall 
become ir regions are protected by photo resist. Then the Nb cap-layer is removed by SF 6 
reactive rf etching, leaving a homogeneous flat NiCu surface, which is then further Ar ion 
etched to partially remove about 1 nm of the F-layer. These areas, in the finished structures, 
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junction 


facets 


a (jum) W(fim) j®\ 


A/cm 2 ) 


|j7|(A/cm 2 ) 


A° (H A} 


(/mi) 


I 


Vc (/AO 


#1 





1 


50 


10 


85 




41 




1.2 


50 


#2 


7T 


1 


50 


10 




35 




65 


0.77 


18 


#3 


0-7T 


2 


25 


10 


85 


35 


41 


65 


1.0 


24 


#4 


0-7T-0 


3 


16.6 


10 


73 


33 


44 


66 


1.0 


23 


#5 


20 X (0-7T-) 


40 


5 


10 


37 


29.5 


62 


70 


3.0 


11.5 


#6 


0-7T disk 


2 


9; 23.5 




4.6 


13.4 


176 


103 


0.29 


6.6 


#7 


0-7T ring 


2 


310 


2.5 


7.3 


2.5 


139 


239 


3.5 


6.8 



TABLE I: Sample parameters at T = 4.5 K: number of facets N, facet length a, junction width W. 
Critical current densities j® and for junctions #3, #4 and #5 were estimated from fits to I C {B). 
A j and Aj refer to the Josephson lengths of the and 7r parts, respectively. They are calculated 
from Eq. @ using the respective critical current densities j® and j*. I is the normalized junction 



length (diameter for #6, circumference for #7), calculated from Eq. (19). The characteristic 



voltage V c = I™ ax /G, where G is the junction conductance, has been inferred by fitting the I—V 



characteristic at maximum critical current /™ ax to the RSJ curve, Eq. (11). For the disk shaped 
0-7T junction the radius r of the inner tt part and total (outer) radius R are quoted instead of a. 
For junction ^7 the facet length a is given by half of the circumference. 

realize the regions, while the non-etched regions are 7r regions. To finish the process, after 
removing the photo resist, a new 40 nm Nb cap-layer is deposited and, after a few more 
photolithographic steps the full structures are completed having a 400 nm thick Nb wiring 
layer, plus contacting leads and insulating layers. The thickness of the F-layer in the devices 
used here is ~ 5 nm and is different for all devices as they come from different places of the 
chip because of a gradient in the F-layer thickness. 

Several sets of 0, tt, 0-tt, O-tt-0 and 20 x (0-7T-) junctions were fabricated in the same 
technological run. The disk shaped and annular samples were fabricated during another 
run. Parameters of the junctions are presented in Tab. [Tj 
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B. Measurement techniques and analysis of LTSEM signal 



For the measurements the samples were mounted on a LTSEM He cryostage and operated 
at a temperature T « 4.5 K. Low pass filters with a cutoff frequency of 12 kHz at 4.2 K, 
mounted directly on the LTSEM cryostage, were used in the current and voltage leads to 
protect the sample from external noise. Magnetic fields of up to 1.2 mT could be applied 
parallel to the substrate plane and thus parallel to the junction barrier layer. We recorded 
I—V characteristics and I C (B). To detect I c we used a voltage criterion V CI ( 0.2 /xV for 
Figs. [3] and [5j 0.5 /xV for Figs. [2] and |3J 1 fjV for all other figures). 

For selected values of magnetic field, LTSEM images were taken by recording the electron- 
beam-induced voltage change 5V(xo,yo) across the junctions (current biased slightly above 
I c ) as a function of the beam-spot coordinates (xo,yo) on the sample surface. The pe- 
riodically blanked electron beam (using /& rs 6.66 kHz, acceleration voltage 10 kV, beam 
current 250 pA), focused onto the sample, causes local heating and thus local changes in 
temperature-dependent parameters like the critical current density j c and conductivity G' 
of the junction. The beam current also adds to the bias current density in the beam spot 
around (xo,yo), but for all measurements reported here the beam current density is several 
orders of magnitude below the typical transport current densities. Thus, this effect will be 
ignored here. The local temperature rise ST depends on the coordinates x, y and z. For our 
SIFS junctions the relevant depth zq is the location of the IF barrier layer, where changes in 
j c and G' affect the I—V characteristics by changing the critical current I c and the junction 
conductance G. We describe the temperature profile within the barrier layer of our junctions 
by a Gaussian distribution 

(x - x ) 2 + (y - y Q ) 2 ) 



5T(x - x , y - y ) = AT exp 



2a 



(3) 



where xq and yo is the position of the center of the e-beam. The LTSEM images pre- 
sented below are reproduced well by simulations using a = 3.5 /zm; this value was used for 
all calculated images shown below and is somewhat larger than for other LTSEM measure- 
ments, presumably due to the relatively thick top Nb layer. Further, from the beam-induced 
changes SI c of the critical current and the measured temperature coefficient dI c /dT, we 
estimate AT m 0.5 K. To a good approximation the beam-induced change of critical cur- 
rent SI c (xo,yo) is proportional to the beam-induced change of the local Josephson current 
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density [35], 5j s (xo,yo) = j c (xa, yo) sin (f)(x a, yo) at I c . To see this we write 
SI c = I c , oa -I c , oS = J(j sm -j s>oS )df 

= J \j c (T + ST) sin 0(T + ST) - j c (T) sin <f>(T)]df. 

(4) 

Here, the subscripts "on" and "off" refer to electron beam switched on and off. The integral 
J (...) c(f has to be taken over the junction area Aj. The local j c depends on the coordinates 
(x, y) via the Gaussian profile of ST(x, y) and possible sample inhomogeneities. In addition, 
j c is different in the and it parts of the junction, with the values of j® and at a given 
temperature. Assuming that the junction is small compared to A j and that a magnetic field 
B is applied in the (x, y) plane, with components B x and B y along x and y, the Josephson 
phase is given by the linear ansatz 

<f>(x, y, <f) ) = 0o + (2tt/$ ) ■ A{B y x + B x y). (5) 

At I c the initial phase </>o is given such that the supercurrent is maximized. For junctions 
having electrode thicknesses larger than the London penetration depth Xl, the effective 
junction thickness is A pa 2Xl pa d! . For our Nb electrodes, using A^ = 90 nm we estimate 
A pa 180 nm. In general, the phase is different in the "on" and "off" states of the beam[43 s 
HI]. When the electron beam disturbs the junction only slightly this difference may be 
neglected and we obtain 

dj c (x,y) 



5L 



sin <p(x,y) ST '(x - x ,yo) 



df. (6) 



dT 

As can be seen in the lower right inset of FigjTJ at least for some of our junctions the 
normalized value 

' dl c 1 \ dj c 1 



(7) 

dT I C J B=Q dT j c 

(assuming a homogeneous j^,jc) i s about constant (pa — 0.2K -1 ) and roughly the same for 
and 7r parts. Note, however, that the latter statement, although valid for the junctions we 
study here, may not always be true. There are cases, e.g. near a temperature driven 0-tt 
transition [TU] where (dj c /dT)/j c of and 7r parts differ strongly in magnitude and perhaps 
even in sign. Assuming a constant value of (dj c /dT)/j c we can further write 

SI c = ^-AT{j s )(x ,y ), (8) 
8 



where we have used the notation 

5T(x - x ,y -y ) 



(js)(x ,yo) 



j c (x,y) sin 0(0 O , x,y) - 



df, (9) 



AT 

where the brackets indicate the convolution of j s with the beam-induced Gaussian temper- 
ature profile Eq. (|3]). When the size of the beam- induced perturbation is small compared to 
the structures to be imaged, we can approximate the Gaussian temperature profile with a 
5-function, and further simplify the above expression to 

djc 1 

5I C « — ~ ATj c (x ,y )sm<f)(4>o,xo,yo)A s , (10) 
al j c 

with spot size A s « 2na 2 , defining an effective area under a 2D Gaussian distribution. 



Eq. ([Top yields 5I C oc j s (x ,yo) = jc(x ,y ) sin0(x o , yo). Thus, by monitoring 5I C , a map of 
j s at I c , including the supercurrent counterflow areas, can be obtained. Note, however, that 
in general the spot size is not small in comparison to the structures imaged. In particular, 
j c sharply changes sign at a 0-tt boundary. Thus, below, we use expression ^ to calculate 
images (j s )(xo,yo) from the simulated supercurrent density j s (x,y) and compare them to 
the LTSEM images. 

To obtain an LTSEM image we do not measure SI C directly (the signal-to-noise ratio would 
be too small for reasonable measurement times which are limited by long term drifts) but 
bias the junctions slightly above its critical current at a given magnetic field and monitor 
the beam-induced voltage change 5V(xo,yo) as a function of the beam position (xo,yo). 
To understand in more detail the corresponding response SV(xo,yo) and the experimental 
requirements to produce a signal proportional to SI c and thus proportional to j s , we first 
note that at the operation temperature the I—V characteristics can be described reasonably 
well by the RSJ model 



V = sgn(/)v / / 2 -/ c ( J B) 2 /G (11) 

for |/| > \I C (B)\ and V = otherwise. Below we will always assume I > and skip sgn (/). 
Examples for a reference junction are shown in Fig. [TJ The I-V characteristics have been 
recorded at B = and at B = 0.33 mT, corresponding to the first side maximum of I C {B). 
Fits to the RSJ curve are shown by lines. Note that different values of G have been chosen 
for the two fits, which, in principle, is unphysical because G should not depend on B. In 
fact, if one fits these I-V characteristics on a large scale one would get equal values of G, 



-60 -40 -20 20 40 60 
V(MV) 



FIG. 1: (Color online) I—V characteristics of a reference junction (#1 in Tab. [j} at T w 4.5 K 
and applied magnetic flux $ = and 3 / 2( ^0> *- e - a t the principal maximum and first side maximum 



of I C (B). Lines correspond to the RSJ curve ((Til, with I c = 428(118) //A and G = 8.6(9.7) ft" 1 



for $ = ( 3 /2$o) Upper left inset shows normalized conductances G/G(4.5K) vs T for a 0, 7r and 
0-7T junction; G(4.5iT) ~ 9.5 fi^ 1 for the and — tt junction, and ~ 9.3 Q — 1 for the 7r junction. 
Lower right inset shows I C (T) at $ = 0, normalized to 4(4.5 K) = 420, 170, and 124 /xA for the 0, 
7r, and 0-7T junction, respectively. All junctions have dimensions of 10 x 50 fim 2 . 



however the region just above I c will not be approximated well, because (11 ) is strictly valid 
only for (3 C = 2ttI c R 2 C/&o = 0. In case of the I—V characteristic for B = we estimate 
that (3 C ~ 0.5 . . . 0.8. Therefore we adopt fits with field-dependent G to reproduce the I-V 
characteristics near J c in the best way. 

When scanning the beam over a junction, which is current-biased slightly above I c , the 
changes 5I C and SG lead to a voltage change 

SV = -%y/ P-I c {BY - h{B)6h -. (12) 
G 2 cy ' Gy/P - I C (B) 2 V 7 

The change in G is related to the temperature rise caused by the electron beam. Similar 
to the case of the critical current, SG(x ,y ) = J df[(dG'/dT)5T(x — Xo,y — y )]. The up- 
per left inset of FigjTJ shows that the relative change (dG/dT)/G = (dG'/dT)/G r is about 
constant for the junctions investigated, with a value of 0.75 K -1 . We, thus, can write 
5G = {dG'/dT)/G> ■ fdf[G'(x,y)5T(x-x ,y-y )) w (dG'/dT)/G' ■ G>(x ,y )ATA s . In 
general, G'(xo,yo) is mainly set by the insulating AI2O3 layer and will not strongly dif- 



fer for the and n parts. Inserting expressions for SI C and SG into (12) we find for the 
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beam-induced voltage change 



SV = ^^ATiFj - F G ) 



where 



and 



dj c 1 A j j c (x ,y )sm(p(xo,y ) 
dT Jc I c {B)yJ[I/I c {B)Y-\ ' 



(13) 



(14) 



(15) 



We emphasize here that these equations rely on the fact that Eq.(ll) provides a good fit 



to the I—V characteristic in the region of interest and should at most be considered as 
semi-quantitative. 

The response due to term Fq is parasitic, if one is interested in spatial variations of 
the supercurrent density. As Fq > 0, it will give a negative and, if spatial variations of 
G'(xq, yo) are small, a basically constant contribution to 5V for the whole junction area (i.e. 
a negative offset). Fi is the response of interest. To make \Fj\ ^> \Fq\ one needs to satisfy 
the condition 



dG' 1 
~dTG' 



dj c 1 



-i 



AjG'(x ,y ) 



h(B) 



G 



1 



[I/I C (B)Y-1 



(16) 



Ajj c (x , yo) sin <p(x , y )G 

When the conductance is about the same for and n parts of the junction, AjG'(x , yo)/G ~ 
1. Further, restricting requirement (16) to coordinates xo,j/o where | sin0(xo,yo)| ~ 1 one 
obtains 

dT j c 



dG' 1 



1 



[I/I C (B)Y-V 



(17) 



dT G'J \dT j c J Ajj c (x , y ) 
with \(dG' /dT)(l/G')(dj c /dT)- 1 j c \ w 3.75 for our junctions (cf. insets of FigJTJ. As we 
will see, when taking images at the maxima of I C (B), at least for Ajj c (xo,yo)/I c (B) w 1, 



Eq. (17) requires the bias current to be less than 10% above I C (B). Note, however, that 



there are cases where Ajj c (xo,yo)/I c (B) is large, e.g., for a homogeneous junction in high 
magnetic field or for a multi-facet junction when the supercurrents of the and tt segments 
almost cancel. In this case the Fq term is not dominant even much above I c . On the 
other hand, to obtain a linear relation between SV and j s (x,y), I should be so far above 
I c that \J '[1 1 1 c(B)} 2 — 1 varies only weakly when the beam is modulated. Typically, this 
requires / to be higher than about 1.05I C (B), leaving only a small window to properly bias 
the junction, i.e. having a response 6V(x ,y ) oc j s (xo,Uo) ■ 
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III. RESULTS 



In this section we discuss I C {B) patterns and LTSEM images of a variety of SIFS junc- 
tions. All data were obtained at T « 4.5 K. For reference, we will start with rectangular 
homogeneous and tc junctions and then turn to rectangular junctions consisting of two, 
three and forty and ir segments. Finally, we will discuss annular and disk shaped 0-tt 
junctions. Sketches of the different geometries are shown as insets in figures |2^a) to[7](a). 

A. Rectangular Junctions 

For all rectangular junctions of length L and width W we use a coordinate system with 
its origin at the center of the junction, so that the barrier (at z — 0) spans from —L/2 to 
+L/2 in x direction and from — W/2 to +W/2 in y direction. 

1. and tt Josephson junctions 

We first discuss results obtained on a junction (#1 in Tab. [T]). Fig. [2] shows I C (B) 
dependence, LTSEM images SV(x, y) and corresponding line scans SV(x) taken at y — 
0. The left hand ordinate of Fig. [2]^a) gives J c in physical units while on the right hand 
ordinate we have J c normalized to I c q = Ajj®. In the graph we compare I C {B) to the 
Fraunhofer dependence, I C {B) = / c (0)| sin(p/(p\, with tp = ttBAL/Qq. In fact, having more 
complex structures in mind, rather than using the analytic expression, we have calculated 
the simulated curve in Fig. |2](a) as 



where (f>(x, y, </>o) is a phase ansatz. Unless stated otherwise, we will assume a linear phase 
ansatz as given by Eq. ([5]). We note here that for junctions containing both and tt segments 
A may differ by some 5 ... 10 % in and 7r regions [1E1 HH] . However, for the sake of simplicity, 
we ignore this effect here. 

For the present junction we have used j c (x,y) = j® = const. The resulting calculated 
I C {B) curve, shown by the dashed line in Fig. [2|a), agrees with the experimental one, con- 
firming the assumed homogeneity of j®. From the value of / c (0) we find j® ~ 85 A/cm 2 and 
Aj ~ 41 fim. Thus, the junction is in the short junction limit with L/Aj ~ 1.2, justifying the 




(18) 
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FIG. 2: (Color online). junction #1 in Tab. |TJ (a) I C (B) with £ || y. Solid (black) curve 



shows experimental data; dashed (red) curve is calculated using Eq. (18). Inset shows the junction 
geometry, (b)-(g) LTSEM images 8V(x. y) |47] taken at bias points indicated in (a), (h)-(m) line 
scans: solid (black) lines 8V(x) are extracted from the corresponding LTSEM images at y = 0; 
dashed (red) lines (j s )(x)/j® are calculated using a ID version of Eq. ([9]). 

use of the linear phase ansatz ([5]). Further, by comparing the abscissas of the experimental 
and simulated curves, one finds that $ = $o corresponds to B w 0.22 mT. From this we 
obtain A »i 200 nm in good agreement with the value of A rs 2A^ ~ 180 nm. Note that due 
to a magnetic field misalignment there will be a slight out-of-plane field component subject 
to flux focusing by large area superconducting films |50j. This leads to an increased value 
of A calculated using the above procedure. 

Fig. |2](b) shows an LTSEM image at B = 0. The corresponding line scan is shown by 
the solid line in Figj2^h). For 5V(x) oc j s (x) = j® one would expect a constant response 
within the junction area. The actual response is somewhat smaller at the junction edges 
than in the interior. Taking the finite LTSEM resolution into account, i.e. calculating 
the convoluted supercurrent density distribution from Eq. (|9j), one obtains the dashed line 
which follows the measured response more closely, although there are still differences that 
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may be caused by the junction, either by a parabolic variation of j® or by a variation in 
conductance G'. To test this we implemented a parabolic variation of along x in the 
calculation of I C (B) / I c (0) and found that the main effect is a slight reduction of the first 
side minima. To still be consistent with the measured I C {B) the variation should be well 
below 10% and is thus most likely not the origin of the 5V variation. To discuss a potential 



G' effect we quantify the SV response using Eq. (13). For the image the bias current was 
set to 1.051c. The function Fq amounts to 0.24 K -1 while for Fj we obtain 0.62 K -1 , i.e. 
changes in conductance contribute by about 1/3 to the total signal. Thus, variations of G' 
in principle could be responsible for the observed variation of SV . However, while we could 
accept a simple gradient of G' along x, the bending in 5V which is symmetric with respect 
to the junction center, is hard to understand. We thus do not have a clear explanation 
for the parabolic shape of 5V(x). To quantify the LTSEM response further, we can look 



at its maximum value AV ~ 0.45 /xV. With I c /G ~ 50 /xV, from Eq. (13) one estimates 
ATA S /Aj pa 0.025 K and from that a beam-induced temperature change AT pa 0.2 K, which 
is somewhat less than 0.5 K estimated from beam-induced I c changes. 

Fig. gc) shows the LTSEM image taken at the first side maximum of I C {B). The field- 
induced sinusoidal variation of 5V(x) can nicely be seen. The corresponding line scan is 
shown by the solid line in Fig. gi) together with (j s )(x), calculated using Eq. Here, a 
potential parabolic-like variation of 5V(x), if present, would be overshadowed by the stronger 
field-induced variation. However, the sinusoidal variation of 5V(x) with an amplitude of 
0.47/xV around an offset value of —0.13 /xV points to beam-induced changes in conductance. 
With the bias current I = 1.1I C (B) we find Fq ~ 0.35/K and Fj pa —1.6/ K, i.e. we expect 
a 20% shift of the sinusoidal supercurrent-induced variation of 5V towards lower voltages, 
roughly in agreement with observation. Further, from the modulation amplitude of 0.47 /xV 
and I c /G = 12 /xV we estimate ATA S /Aj ~ 0.025 K in agreement with the estimates for the 
zero field case. 

Finally, Figs, pfd)-(g) show LTSEM images and Figs. j2j(j)-(rn) corresponding line scans 
for higher order maxima in I C (B). In all cases, the field-induced modulation of SV(x) can be 
seen clearly, and simulated curves for (j s )(x), calculated using Eq. (|9]), are in good agreement 
with measurements. 

We found similar results also for other reference junctions, including 7r ones. In the latter 
case, typical values at T pa 4.5 K of the critical current densities are j* ~ 30 A/cm 2 (see e.g. 
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FIG. 3: (Color online). 0-tt junction ^3 in Tab. |T} (a) I C (B) with B \\ y. Solid (black) curve 



shows experimental data; dashed (red) curve is calculated using Eq. (18). Inset shows the junction 
geometry. (b),(c) LTSEM images 5V(x, y) (37] taken at bias points indicated in (a). (d),(e) corre- 
sponding images (j s )(x,y)/j® calculated using Eq. Q. (f),(g) line scans: solid (black) lines 5V(x) 
are extracted from the corresponding LTSEM images at y = 0; dashed (red) lines (j s }(x)/j® are 
calculated using a ID version of Eq. (|9j). 



#2 in Tab. |Ij). This value is not large, but it is almost an order of magnitude higher than 
what has been previously reported for SIFS tt junctions [15]. 



2. 0-tt Josephson junction 

Now we discuss data for a 0-tt junction (#3 in Tab. [T]) presented in Fig. [3] The simulated 
I C (B) curve in Fig. [3]^a) fits the experimentally measured dependence in the best way for 
Jc /Jc = ~ 0.42. The right hand axis is normalized to I c q = j®Aj. From the measured value 
of / c (0) and the junction area Aj we find j° = 85 A/cm 2 and = —35 A/cm 2 . For a 0-tt 
junction, Aj can only be defined in and tt parts separately, but not for the junction as a 
whole. However, one can find a normalized junction length as 
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where L and L n are the total lengths of and it parts and Aj and A} are the Josephson 
lengths in the and 7r parts, respectively. With this definition we calculate I ~ 1, showing 
that the junction is again in the short limit. For A we obtain a reasonable value of 200 nm. 
Further note that the measured I C {B) is slightly asymmetric, i.e. the main maximum at 
negative field is slightly lower than at positive field. This effect, which is not reproduced by 
the simulated curve, is due to the finite magnetization of the F-layer which, in addition, is 
different in the and 7T parts. This effect is addressed elsewhere [48J. 

For the 0-tt junction, at B=0 the supercurrents of the two halves should have opposite 
sign. The part giving the smaller contribution to J c should show inverse flow of supercurrent 
with respect to the applied bias current, i.e., the 7r part in our case. This can be seen 
nicely in Fig. |3|b) showing an LTSEM SV(x,y) image at zero field. The 7r part is on the 
left hand side. For comparison, Fig. [3^d) shows a (j s )(x,y)/ j® image of the supercurrent 
density distribution, calculated using Eq. For better comparison, Fig. [3]T) shows a 
measured and a calculated line scan. The left ordinate is shifted by 0.47/iV relative to the 
origin of the right ordinate to match the simulated and experimental curves. This shift is 
required to account for the beam-induced conductance change. More quantitatively, with 
J/J c (0) ~ 1.06, Iq/G = 13.5/iV and assuming that G' is the same for and 7r parts, we 
estimate Fq ~ 0.3 K _1 . For the tc part we estimate F/ w — 0.8K -1 , while for the part 
we obtain Fj m 1.9K -1 . The peak-to-peak voltage modulation in the LTSEM image is 
1.65 /iV. From these numbers we estimate (A s /Aj)AT ~ 0.045 K, or AT « 0.3 K, which is 
reasonable. For the conductance-induced shift we obtain a value of about —0.2 /xV, which 
is about a factor of 2 less than expected from the measurement, but still within the error 
bars. 

The LTSEM image 8V(x, y) shown in Fig. |3](c) has been taken at the main maximum 
of I C (B). Here, both parts of the junction give a positive response. The measurement is in 
good agreement with expectations, as can be seen in the calculated image (j s )(x,y)/ j® in 
Fig. [3^e) and by comparing the line scans 5V(x) and (j s )(%)/jc shown in Fig. [3^g). Note 
that the "offset problem" seems to be less severe here. Indeed, with I c /G = 24 /iV and 
I/I c = 1.019 we obtain F G « 0.15 K" 1 and Fj « 2.1 Kr 1 for the part and Fj « 0.85 KT 1 
for the 7r part. The supercurrent term thus clearly dominates. 
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FIG. 4: (Color online). O-vr-0 junction #4 in Tab. |T} (a) 4(5) with B \\ y. Solid (black) curve 



shows experimental data; dashed (red) curve is calculated using Eq. (18). Inset shows the junction 
geometry. (b),(c) LTSEM 5V(x,y) images [4"T] taken at bias points indicated in (a). (d),(e) corre- 
sponding images (j s )(x,y)/ j® calculated using Eq. Q. (f),(g) line scans: solid (black) lines 5V{x) 
are extracted from the corresponding LTSEM images at y = 0; dashed (red) lines {js)( x )/jc are 
calculated using a ID version of Eq. (|9|). 



3. 0-7T-0 Josephson junction 

Next we discuss data for a 0-7T-0 junction ($4 in Tab. [T]) presented in Fig. |1J The best fit 
to I C (B) was obtained for = 73 A/cm 2 and = —33 A/cm 2 . From here we obtain(/ w 1). 
We are thus again in the short junction limit. Further, we obtain A w 200 nm. 

LTSEM images, taken at, respectively, the central maximum and the main maximum at 
positive fields, are shown in Figs. Efb) and (c). Figs. |4](d) and (e) are simulated images, 
and Figs. |i](f) and (g) show the corresponding line scans. For this junction, the simulated 
curves, taking only modulations due to j s into account, agree well with the data. For 
Fig. gb), with I/I e = 1.044 and IJG = 17.5 we find F G w 0.22 K" 1 and, for the j s 
maximum in the part, Fj w 1.55 K -1 . For the j s maximum in the 7r part we obtain 
Ff ~ 0.7 K -1 . The offset is thus not very large. From the peak-to-peak modulation of 
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1.35 /iV we estimate (A s /Aj)AT pa 0.035 K and, thus, a reasonable value AT pa 0.23 K. 
Taking this value, we estimate the offset voltage to about 0.1 /xV. For the measurement at 
the main maximum with I/I c = 1.04 we obtain IJG = 23 /xV, F G pa 0.21 K -1 , F? pa 1.29 K" 1 
and Fj pa 0.58 K -1 . Using (A s /Aj)AT = 0.035 K we expect an offset in 5V of -0.17/xV 
and a maximum supercurrent response of 0.85 /xV in the parts, and 0.3 /xV in the central 
7r part. The measured numbers are 0.65 /xV and 0.35 /xV, respectively. 



4- 20 x (0-7T-) Josephson junction 

Having seen that well behaving 0-7T-0 junctions can be fabricated one may consider mul- 
tisegment structures where many 0-7r segments are joined. The main purpose here is to 
check the complexity and reliability of the structures that can be fabricated already now. 
Moreover, as already mentioned in the introduction, multi-segment iV x (0-7T-) Josephson 
junctions are promising for the realization of a <p junction. The structure we study here has 
twenty 0-n segments (#5 in Tab.JT]). In Fig. |5|a) we compare the measured I C (B) dependence 



(solid line) with the one calculated (dashed line) using Eq. (18) with a linear phase ansatz 
(pi). However, on both sides of each main peak we see quite substantial deviations of the 
calculated curve from the experimental one. In particular, the series of J c maxima following 



the main peak are much higher in experiment than in simulations based on Eqs. (18) and 
(pp. It is interesting that such a shape of I C (B) was also measured for d-w&ve / s-w&ve zigzag 
shaped ramp junctions[32| |33| 134"] . 

To understand the origin of such deviations, we have tested numerically a variety of local 
inhomogeneities j c (x) in the different facets, ranging from random scattering to gradients 
and parabolic profiles, always using the linear phase ansatz (|5]). None of them, and also 
no variations in effective junction thickness A(x) were able to qualitatively reproduce the 
I C (B) features described above. Finally, it turned out that the quantity to be modified is 
the phase ansatz, i.e., the field becomes non-uniform. Adding a cubic term, which accounts 
for a small phase bending, we have (assuming B || y) 

3" 



B V AL 

(p(x, y, O ) = 0o + 2tt— — 



2x (2x 



(20) 



Calculating I C (B) using Eq. (18) with from Eq. (20), we were able to reproduce the above 



mentioned features of the experimental I C {B) dependence, as shown by the dotted line in 
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FIG. 5: (Color online). 20 x (0-tt-) junction #5 in Tab. ^ (a) I C (B) with B \\ y. Solid (black) 



curve shows experimental data; dashed (red) curve is calculated using Eq. (18) with linear phase 
ansatz (|5]); dotted (blue) curve is calculated using Eq. (18) with cubic phase ansatz (20). Inset 
shows the junction geometry. Only 10 0-7T- segments are drawn for simplicity, (b) LTSEM images 
SV(x, y) |47j taken at the bias points indicated in (a) by numbers 1 to 8. (c) corresponding images 



(j s )(x,y)/jc calculated using Eq. ^ and the cubic phase ansatz (20). (d) line scans: solid (black) 
lines 5V(x) are extracted from the corresponding LTSEM images at y = 0; dashed (red) lines 
(js)(x)/jc are calculated using Eq. 

Fig. [5^ a). Here we used 03 = —0.065, i.e. a rather small correction to the linear phase. In 
spite of this, for the relatively high magnetic fields around the main maxima of I C {B), this 
term adds up to an additional phase ~ 1 and becomes important — the contribution to the 



integral in Eq. (18) changes essentially close to the junction ends. Note that a homogeneous 
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junction or a junction consisting of only a few and ir segments could not sense that, since 
at the high fields, where the bending of the phase reaches values of ~ 1 at the junction 
edges, I c is already suppressed to almost zero. 

As we will show in a separate publication [49] the origin of the nonlinear contribution in 



Eq.(20) is a parasitic magnetic field component perpendicular to the junction plane, which 
appears due to a misalignment ~ 1° between the (x, y) plane and the applied magnetic 
field. This perpendicular component causes screening currents that result in a non-uniform 
(constant+parabolic) field focused inside the junction and pointing in y direction. Similar 
effects can also be present in non-local planar junctions|51j. but we are far from this limit. 

By comparing the nonlinear-phase simulation to the measured I C {B) we infer j® = 
37A/cm 2 , = — 29.5A/cm 2 and I ~ 3. The junction is thus still in the short limit. 
We further obtain A « 350 nm, which is higher than the value we obtained for the other 
rectangular structures, but consistent with the fact that we have a focused out-of-plane field 
component. 

Fig. [5^b) shows a series of LTSEM images. Image 1 is taken at B = 0, image 2 at the 
main maximum and images 3 to 8 at the subsequent maxima. For image 1 one can nicely 
see the modulation induced by the 40 facets, although negative signals are not reached any 
more. This is due to the small facet size of 5 /xm which is on the LTSEM resolution limit. 
At the main maximum the signal is strong and positive, with a slight long-range modulation 
but no evidence of modulations due to the individual facets any more. At the higher maxima 
(images 3 to 8) additional minima appear in 6V(x, y). Fig. |5]^c) shows the corresponding 
images calculated using the cubic phase ansatz, and Fig. ^d) shows the corresponding line 
scans, comparing the measured 8V(x) (solid lines) with the calculated {j s ){x) (dotted lines). 
As can be seen, the agreement is excellent, except for the line scan taken at the I c maximum. 
Here, the measured response is strongly weakened towards the junction edges in contrast to 
the calculated modulation of j s . For this bias, with / = 1.029/ c we estimate Fq ~ 0.18 K -1 
and Fi « 1.5 K . It is thus not very likely that the discrepancy is caused by a spatially 
varying conductance. On the other hand, from the well behaved LTSEM images at zero field 
we can rule out a long range variation of j% and j* as well. A possible origin of this behavior 
may be a non-uniform field focusing that results in a phase ansatz <p(x, y, <po), which is more 



complicated than the cubic one of Eq. (20). However, we have to admit that we did not 
succeed in finding a proper dependence. 
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FIG. 6: (Color online). Disk shaped 0-7r junction #6 in Tab. [Tj (a) I C {B) with -B || y. Solid (black) 



curve shows experimental data; dashed (red) curve is calculated using Eq. (18); dotted (blue) curve 



is calculated using Eq. (21 ). Inset shows the junction geometry, (b)-(d) LTSEM images 8V(x, v) |47j 
taken at bias points indicated in (a), (e)-(g) corresponding images (j s )(x,y)/j® calculated using 
Eq. ([9]). (h)-(j) line scans: solid (black) lines SV(x) are extracted from the corresponding LTSEM 
images at y = 0; dashed (red) lines (j s )(x)/jc are calculated using Eq. pi). 



We have measured several 20 x (0-7T-) junctions. All behaved similar to the one discussed 
here, including the shape of I C {B) with a well developed set of maxima following the main 
peak and also with respect to LTSEM images. Thus, the present SIFS technology is fully 
able to deliberately produce quite complicated multi-facet 0-7r junctions. 
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B. Disk Shaped and Annular Junctions 



1. Disk shaped Jo sephson junction 

The SIFS technology offers the possibility to create a more complex 0-ir boundary than a 
linear one. An intriguing option is to close this boundary in a loop. The disk shaped junction 
#6 in Tab. [I] is of this type. Here, we use a coordinate system with its origin at the center 
of the disk, see the inset of Fig. [6|a). The I C (B) dependence, shown in Fig. [6|a), exhibits a 
central maximum at B = where the critical currents of the and the tt part subtract, as 



well as prominent side maxima. By fitting the curve calculated using Eq. (18) (dashed line) 
to the experimental curve (solid line), we obtain j% = 4.6 A/cm 2 and j* = — 13.4A/cm 2 as 
optimal fitting parameters. Referring to 2R as the junction length we obtain I « 0.29, i.e. 
again the junction is in the short limit. Fitting the horizontal axis using the length 2R we 
obtain A ps 200 nm. 

For this sample, I c /G ~ 6.8 /xV (at zero field) is rather low. As a consequence the 
detectability of I C {B) at low values of the critical current is resolution limited. We used 
a voltage criterion V cr = 1 /xV to measure the "critical current" , yielding a parasitic J c 
background of I r = V CI G ~ 6 /xA. When comparing simulation with experiment the value 
of I r should be added (in quadrature) to the calculated critical current If m to obtain the 
"visible critical current "7™, which should be compared with the experimental one I* xp , i.e., 



(i^y + (i r y. (21) 

One can see in Fig. |6|a) that the calculated curve including I r (dotted line) is in good 
agreement with the experimental data. 

Fig. |6^b) shows an LTSEM image 5V(x,y) taken at the central maximum of I C (B). 
Fig. [6|e) shows the corresponding simulation of (j s )(x,y)/j^ and Fig. [6^h) contains cor- 
responding experimental and calculated line scans. The LTSEM data and the simulation 
results agree well, showing that the supercurrent in the central tc region flows against the 
bias current. Figs. [6^c),(f),(i) show the results for an applied magnetic field corresponding 
to the first side maximum of the I C {B) curve. Here, the field-induced sinusoidal variation 
of the supercurrent is superimposed with the disk shaped 0-tt variation. The supercurrents 
in the ir region as well as in a major part of the region flow in the direction of the bias 
current, maximizing I c . For completeness, in Figs. |6](d),(g),(j) we also show corresponding 
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plots taken at the second side maximum of I C (B). Here, the magnetic field induces about 7 
half oscillations of the supercurrent density along x. Similar to the previous cases, experi- 
mental and calculated plots agree well. For the central maximum with I/I c = 1.09 we find 
Fq ~ 0.3 K -1 and F® pa 2.5 K _1 , FJ pa 7.2 K . Thus, the offset due to conductance changes 
is minor in this case. The same holds for the other bias points. The main reason is that the 
factor j c (xo,yo)Aj/ I C (B) entering Fi is large (e.g. about 7 for the 7r part at B = 0). 

2. Annular Josephson junction 

The last structure we want to discuss in this paper is an annular 0-tt junction (#7 in 
Tab. [TJ see the sketch in Fig. [7|. Half of the ring is a region and the other half is a 7r 
region. One thus obtains an annular junction with two 0-tc boundaries. If the junction were 
long in units of Aj it would be a highly interesting object to study (semi)fluxon physics, 
similar to the case of Nb junctions equipped with injectors[52j [53] . For this junction we use 
a coordinate system with its origin in the center of the ring, and the steps in the F-layer 
are located on the y axis. Fig. (7^a) shows I C (B) of this structure, with B \\ y. The critical 
current is always above 10 fiA. This offset is in fact real and reproduced by the simulated 
I C (B) which is for I r = (the actual value I r ~ 8 fiA only slightly lifts the I C (B) minima). 
From the fit we obtain a ratio j^/j® = —0.35. Taking into account that Aj pa 1550 fim 2 , 
we get j® pa 7.3A/cm 2 and ~ — 2.5A/cm 2 and, referring to the circumrefence as the 
junction length, I pa 3.5. Thus, we are still in the short junction limit. Further, we obtain 
A pa 150 nm, which is somewhat lower than for the other junctions, but still reasonable. 

Figures [7j^b)-(e) show LTSEM images taken at various values of B as labeled in the I C (B) 
pattern shown in Fig. [Wa) . As shown in Fig. ^b) for B = 0, i. e. at the central local minimum 
in I C (B), a counterflow in the tt part (left half) can be seen. At the main I c maximum the 
supercurrents in both the and the 7r region flow in the direction of bias current [Fig. [7](c)]. 
Images (d) and (e), taken at the subsequent I C (B) maxima, look more complicated, showing 
several regions of counterflow. In all cases, however, the LTSEM images are well reproduced 
by simulations, as can be seen in Figs.[7](f)-(i) and the corresponding linescans, see Figs.^j)- 
(m). The linescans, taken along the junction circumference, start at the upper 0-7r boundary 
and continue clockwise. 

For this annular junction we have also rotated the magnetic field by about 70° towards 
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FIG. 7: (Color online). Annular 0-7T junction #7 in Tab. [T] with B \\ y: (a) I C {B) pattern; solid 



(black) curve shows experimental data; dashed (red) curve is calculated using Eq. (18). Inset shows 
the junction geometry, (b)-(e) LTSEM images jJTj taken at bias points indicated in (a), (f)-(i) 
corresponding images calculated using Eq. (|9]). (g)-(m) line scans: solid (black) lines 5V(x) are 
extracted from the corresponding LTSEM images; dashed (red) lines (j s ){x)/jc are calculated using 
Eq. ([9]) with the curvilinear coordinate x instead of x which runs along the junction circumference. 



the x direction. The corresponding data are shown in Fig. [8j For this field orientation I C (B) 
strongly differs from the case B \\ y, cf., Fig. [8](a), but can be reproduced by simulations, 
using the same j° and j* as in Fig. [Tj Furthermore, simulations show that if the field 
is rotated further towards the x axis, the height of the side maxima in I C {B) decreases, 
reaching only half of their height of the 70° case when the field is parallel to the x axis 
and the I c minima reach zero. Thus, the annular 0-7r junction reacts very sensitive to field 
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FIG. 8: (Color online). Annular 0-n junction #7 in Tab. [I] with B applied in the (x,y) plane 
under 70° from the y direction: (a) I C {B) pattern; solid (black) curve shows experimental data; 
dashed (red) curve is calculated using Eq. ( |18[ ). Inset shows an SEM image of the junction, (b) 
and (c) LTSEM images [TF] taken at bias points indicated in (a), (d) and (e) corresponding images 
calculated using Eq. (|9]). (f) and (g) line scans: solid (black) lines SV(x) are extracted from the 
corresponding LTSEM images; dashed (red) lines (j s )(x)/jc are calculated using Eq. ^ for x. 

misalignments relative to the x axis, similar to the case of the 20 x (0-7T-) junction where 
out-of-plane field components strongly altered I C (B). For completeness, Fig. |8](b)-(g) also 
shows LTSEM images taken at the selected bias points labeled in Fig. [8](a) and compare 
them with simulation. The agreement is again very good. 

IV. CONCLUSION 

We have studied a variety of SIFS Josephson junction geometries: rectangular 0, it, 0-ir, 
0-7T-0 and 20 x (0-7T-) junctions, disk-shaped 0-tt junction, where the 0-ir boundary forms 
a ring, and an annular junction with two 0-tt boundaries. Using LTSEM we were able to 
image the supercurrent flow in these junctions and we demonstrate that and 7r parts work 
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as predicted having j® > and < 0. Within each or tc part, according to both LTSEM 
images and I C (B), the critical current density is rather homogeneous. Particularly, within 
our experimental resolution of a few /iin, we saw no inhomogeneities that might have been 
caused by an inhomogeneous magnetization of the F-layer. This implies that ferromagnetic 
domains, although probably present, must have a size well below 3 /mi. 

These results demonstrate the capabilities of the state-of-the-art SIFS technology. Ar- 
rangements like the ring-shaped 0-tc boundary are impossible to realize using other known 
0-7r junction technologies [3UJ E2J EHl Even intersecting 0-tt boundaries seem to be fea- 
sible, e.g., by arranging and it segments in a checkerboard pattern. 

For the 7r regions we demonstrated a record value of ~ 35 A/cm 2 at T « 4.5 K, which 
is an order of magnitude higher than the values previously reported for SIFS junctions with 
a NiCu F-layer [161 EI]- Still, to obtain reasonable values of Xj < 20 /mi, should be 
increased by at least one order of magnitude to reach ~ lkA/cm 2 . Then the 0-7r junctions 
can be made long enough (in units of Aj) to study the dynamics of semifluxons pinned 
at the 0-7r boundaries. In this case semifTuxon shapes, not realizable with other types of 
junctions, are possible, e.g., closed loops, intersecting vortices, etc.. Another issue inherent 
to the present SIFS technology is that the critical current densities j° and in the and ir 
parts are not identical in general. In many cases this does not matter, e.g., when one works 
with semifluxons in a long junction. If = is required, the difference in j® and will 
lead to a low yield of the circuit and one is perhaps restricted to operate the device in a 
narrow temperature interval, where j° and are closer to each other. 

Even with the present restrictions, quite complex geometries like the 20 x (0-7T-) junctions 
have been realized. Those SIFS multifacet junctions already showed interesting features, 
like their high sensitivity to nonuniform magnetic fields, and they will be usable for many 
fundamental studies, e.g. on the way of realizing if junctions. 
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